Gapless Fermi Surfaces in Superconducting CeCoIns? 
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According to [M.A. Tanatar et al, Phys. Rev. Lett. 95, 067002 (2005)], in a multi-band d- 
wave superconductor CeCoIns électrons remain partially uncondensed. Interactions must induce 
superconducting order on ail Fermi surfaces. We calculate spécifie heat and thermal conductivity 
in a two band model in présence of defects. Superconductivity originates on one Fermi surface, 
inducing a smaller gap on the other. Impurities diminish the induced gap and increase the density 
of states, restoring rapidly the Wiedemann- Franz law for this Fermi surface. Our calculations are 
in agreement with experiment. 
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The heavy fermion "115" materials CeMIn5 
(M=Ir,Co,Rh) are currently under extensive expér- 
imental scrutirw. The compounds reveai non-Fermi 
liquid features[ï|, the "two-fluid" , or, alternatively, 
two energy scales behavior in their thermodynamics 
and magnéto-transport 0. Thèse materials are d-wave 
superconductors (SC), and their properties in high 
magnetic fields testify in favor of an important rôle 
played by the Pauli spin limitation and, more generally, 
the involvement of magnetic mechanismJ3|. 

More récent experiments on thermal conductivity 
and spécifie heat in the heavy fermion superconductor 
CeCoInsJJl have found that the variation of thermal 
conductivity with increased concentration of impurities 
(Cei_a;La2;CoIn5) contradicts the expected behavior of 
nodal quasiparticles in a d-wave superconductor 1^. 
According to Ref.Q, thèse results favor multi-band su- 
perconductivity, with a group of light quasiparticles that 
remains uncondensed and coexists with superconductiv- 
ity on other Fermi surfaces (FS), even in the absence of 
defects {x — 0). It was suggestedQ that such a scé- 
nario could explain many other expérimental anomalies 
observed in this material, such as the saturation of the 
pénétration depth at low temperatureslâ] , the unexpect- 
edly large différence in thermal conductivity^ between 
longitudinal and transverse field orientations in the plane 
with respect to the heat flow, and the apparent contra- 
diction in the angular field dependence of thermal con- 
ductivity and spécifie heatp^. 

It seems difficult to justify the idea of uncondensed 
électrons on separate FS in the SC stateT.lJ. There are 
no fundamental reasons why the interactions between 
électrons on any two given Fermi surfaces in a multi- 
band métal would become identically equal to zéro. In 
other words, if the superconducting (SC) pairing takes 
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place on one of the FS, even an arbitrary weak inter- 
action would inevitably induce a SC order parameter 
on the second FS[Î2. At the same time, it is well- 
known that ordinary impurities may lead to gapless su- 
perconductivity in multi-band superconductors. One ex- 
ample is MgB2, where at least two différent Fermi sur- 
faces participate in s- wave superconductivity|l^. An- 
other example of a multi-band case is given by surface 
superconductivity 13, 'ÏE'l, where a 2D Fermi surface be- 
comes split by the Rashba spin-orbit interaction 16, 17]. 

The "115" CeMIng materials (where M = Rh,Ir,Co) 
are multi-band compounds. In particular, in CeCoIns 
quasi-2D Fermi surfaces coexist with small 3D onesflSl 
hj. CeCoIns is a d-wave superconductor [lÔj . 

The main assertion made in Réf. 0] is that while in 
the d-wave Cei-ajLa^jCoIns the thermal conductivity, k, 
must obey the so-called "universal" law[^|^, experimen- 
tally it does not, and gapless carriers are needed to ex- 
plain the observed behavior. 

According to the "universal" law, the zero-T limit 
value of {k/T)t=o should not dépend on the concentra- 
tion of defects, x. The argumentation given in Ref.0] 
has foUowed the ideas of the kinetic theory. It is well 
known that in the absence of defects the density of states 
(DOS) at the d-node is zéro at zéro energy and increases 
linearly with energy [23, |^. Defects would introduce a 
finite DOS. If the latter were proportional to x, spécifie 
heat would be proportional to x as well. With the num- 
ber of scattering centers given by the very same x, the 
concentration of defects drops out of the new relation 
that connects the thermal conductivity and resistivity 
instead of the familiar Wiedemann- Franz law. 

Actually, the "universal" law has been rigorously 
derivedla, "o] only in the so-called "clean limit" . By this 
one means that the concentration of defects, x, is very 
small, so that one can neglect the changes in Te and in 
the value of the gap order parameter with x. 

Another provision for the validity of the arguments 
given in Ref.Q is that, at the same time, x should be 
not too small, otherwise one would encounter the région 
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of a well-known singularity in DOS that appears due to 
the présence of impurity scattering near x = 0[2^. where 
the spécifie heat (or DOS) cannot be linear in x. The 
range of applicability of this second provision is not well 
defined, and, as we show below, dépends on the strength 
of impurity centers. 

The arguments given in Ref.0| may look plausible in 
the framework of the kinetic theory, and may even be 
fulfilled for the Zn-doped cuprates[3,la,l23i ^'^s certainly 
not applicable for the samples of Cei-j^La^^ Coins studied 
in Réf. 3. In fact, the latter are not in the "clean limit", 
as foUows from Figs 1,2 of Réf. 3, where at a; = 0.02 
the value of Te is already significantly reduced 0, and, 
hence, r^^ ~ Te ^ A. For ail the above reasons the 
data presented in Fig.3 of Réf. Q for the thermal con- 
ductivity k(T,x)/T in the alloy Cei-a^La^CoIns cannot 
be analyzed in terms of the "universal" law. Note, how- 
ever, that at the lowest x the data j] reveals the prés- 
ence of a thermal conductivity component that rapidly 
decreases with the increased x, the feature reminiscent of 
the x-dependence for normal phase thermal conductivity, 
where the Wiedemann- Franz law applies. 

To accommodate this last fact, we postulate the exis- 
tence of the second FS (FS2) with a small induced SC 
gap, which is therefore very sensitive to scattering on de- 
fects. This FS, we suggest, becomes practically gapless 
already in the présence of a small amount of defects, with- 
out changing significantly SC parameters on the other, 
leading FS (FSl). FSl makes its own contributions to 
the spécifie heat and thermal conductivity. We argue 
that the sum of contributions from the two FS explains 
the characteristic z-dependence for dataPf without in- 
voking an unphysical idea of isolated and uncondensed 
électrons on one of the FS. 

In what follows we perform calculations of the spé- 
cifie heat and thermal conductivity for two FS separately: 
first, for FS2, where one may negleet the variation of the 
order parameter with the concentration of defects, and, 
second, for FSl, where such variations must be taken into 
account. The model, as it is shown below, accounts well 
for the main feature in the expérimental data 4] , at least 
qualitatively. 

We use the standard weak-couplin g BC S-like approach 
to unconventional superconductivit vl2ll l24l. generalized 
to the muhi-band situation^ [H |2j|2||. The Fermi 
liquid quasiparticles in bands i and j interact via a weak 



short-range interaction, 
be written as: 



Uij{r,r'). The Hamiltonian can 
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Hère cr, tr' are spin indices, and the sum over ail indices 
is written in momentum space, k, fc', q. The scattering 
potential for impurities has the form: 

= E E E e'^''-''^''-Vi,{k, q)al^ai,„ (2) 

where l and j are band indices and Ra are the defect po- 
sitions. In turn, the inverse Gor'kov matrix defining 



the normal (G) and anomalous {F, F^) Green's function 
is: 
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where, as usualj29l|. iûij,n = i(^n — '^ji'^n)- In the Born 
approximation: 
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is the concentration of impurities, while the self- 
energy for the anomalous Green's function Fj{ioJnk) van- 
ishes for isotropic scattering due to momentum intégra- 
tion, and \u,i\^ ^ J ^ J fL |y^,(fef , fc^p. The Gor'kov 
équations are then solved, along with the gap équation. 
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We find: 
where 
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The gênerai analysis of the above équations when 
ail couplings are of the same order is somewhat 
cumbersome fï^ |2^ . Neither is it necessary, since Te dé- 
pends exponentially on the coupling constant, and su- 
perconductivity will first émerge on the Fermi surface 
with the largest quasiparticle attraction potential (de- 
noted FSl below). Even if \Uii\ > IC/22I, |t/i2|, the su- 
perconducting gap induced by the intra-band coupling 
on other Fermi surfaces (FS2, FS3, etc.), although being 
small, has the smallness only linear (not exponential!) in 
Ui2. Clearly, for singlet unconventional superconductors 
in présence of pair-breakers the smallcr gap will get com- 
pletely "erased" first, resulting in a new type of super- 
conducting state, in which almost "normal" (t^^ ^ A2) 
quasiparticles on the FS2 coexist with "gapped" (i.e., 
having the d-wave nodes) quasiparticles on the leading 
band. 

We consider first the case when the amount of defects, 
while being too small to significantly change Te and SC 
parameters on FSl, is sufficient to change thermodynam- 
ics and transport attributable to FS2. 

The theoretical formalism is well-known (see, for in- 
stance, Réf. !3). However, numerical results are différ- 
ent for FS2 (compare, e.g., with Ref.js^]), because A for 
FS2 (= A2) is taken as an x-independent constant. 

For a better readability, we provide a brief summary of 
the theoretical results here. In Eq.Q the analytical con- 
tinuation of z(D„ from the upper energy half-plane defines 
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the variable t{E) (in the Born approximation) as: 

r2iT 



t 



Afc — Acos(2(^). 
The density of states is: 

N{E) = 2NoTlm[t{E)l 
and for the spécifie heat one has: 
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Hère A'o is the normal DOS. From gênerai expressions 
for the thermal conductivity, k(T, ce) la lai, we need only 
the low température limit|y| (see Ref.|2J|): 
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If in Ea. H12|l 7 <C A, one obtains the notorious "uni- 
versal" limit,_^.e., n independent of the concentration of 
impurities[~ 



The formai définition of 7 is given by: 
j = Im[tiE ^ 0)], 



(13) 
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however, its explicit expression strongly dépends on the 
assumption about the strength of the scattering potential 
for the defects|3lj. At small x, in the Born approxima- 
tion, 
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In the opposite limit of strong scattering (the unitary 
limitlâïl): 



where 
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The results of our calculations of the spécifie heat and 
thermal conductivity contributions from FS2 are pre- 
sented in FigH for thèse two extrême cases. One sees 
that in the Born approximation a flat région (the "uni- 
versal" law) in k/T appears at small x; in the unitary 
limit for impurity scattering this interval is practically 
absent, although, as it was explained above, A2 was kept 
constant. At larger x k/T decreases, merging rapidly 
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FIG. 1: Contribution of FS2 to the linear coefficient for the 
heat capacity 7(T)/7Jv(r) (left axis) and the ratio of thermal 
conductivity k{T)/T to its universal limit for FS2 (right axis) 
at small températures T ^ as a function of l/2rA in the 
Born approximation, or F/ A in the unitary limit. 
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FIG. 2: Contribution of FSl to the linear coefficient for the 
heat capacity 7(r)/7jv(r) (left axis; blue-Born approxima- 
tion, green-unitary limit) and the ratio of thermal conductiv- 
ity K,{T)/T to its universal limit for FSl (right axis; black- 
Born approximation, red-unitary limit ) at small tempéra- 
tures T — > as a function of l/2rA in the Born approxima- 
tion, or F/ A in the unitary limit. The inset shows the de- 
pendence of Te (black) and the order parameter in the Born 
approximation (red) and the unitary limit (blue). 



with the 1/x behavior expected from the Wiedemann- 
Franz law in the normal phase. No interval is seen for 
the linear in x behavior of the Sommerfeld coefhcient in 
the spécifie heat. 

The same charactcristics have been calculated for FSl 
and are given in Fig[21 Unlike for FS2, now both A and 
Te decrease with shown in the inset (the dependen- 
cies practically coïncide with the results for the unitary 
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limit from Réf. f3ÔI| ) . Note a considérable variation oî k/T 
both in the Born approximation and in the unitary limit 
for the scattering amplitude between x = and x = Xcn 
whcre Te — and SC is destroyed. A linear in x behavior 
in the spécifie heat is seen in some a;-interval only in the 
unitary limit. Above Xcr DOS rcaches its normal value, 
independent of x, and k/T begins to decrease as l/x, 
according to the Wiedemann-Franz law, leading to the 
corresponding "bump" at X — Xcr- Note that one secs 
such a "bump" in expérimental data shown on Fig.3 of 
Ref.i. 

To summarize, the model of two FS with différent val- 
ues of SC order parameter provides an excellent qualita- 



tive interprétation of the remarkable findings in Ref.|j]. 
Indeed, we have shown that at smaller x the x-dependent 
contribution to k/T in R éf. M can be attributed to a FS2 
(such as the e- sheetfï^ IlOf) "stripped" of a SC gap by 
strong scattering. Linear in x dependence of the Som- 
merfeld coefficient in the spécifie heat at T = seen in 
Ref.^éj may corne from the FS 1*331. 
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